The fractional volatility model: An 
agent-based interpretation 

R. Vilela Mendes* 
Complexo Interdisciplinar, 
Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal 



Abstract 

Based on criteria of mathematical simplicity and consistency with em- 
pirical market data, a model with volatility driven by fractional noise has 
been constructed which provides a fairly accurate mathematical parametriza- 
tion of the data. Here, some features of the model are discussed and, using 
agent-based models, one tries to find which agent strategies and (or) prop- 
erties of the financial institutions might be responsible for the features of 
the fractional volatility model. 

Keywords: Fractional volatility, Statistics of returns, Option pricing, Agent- 
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1 Introduction 

Classical Mathematical Finance has, for a long time, been based on the as- 
sumption that the price process of market securities may be approximated by 
geometric Brownian motion 

dSt = fiStdt + aStdB{t) (1) 

In liquid markets the autocorrelation of price changes decays to negligible values 
in a few minutes, consistent with the absence of long term statistical arbitrage. 
Geometric Brownian motion models this lack of memory, although it does not 
reproduce the empirical leptokurtosis. On the other hand, nonlinear functions 
of the returns exhibit significant positive autocorrelation. For example, there is 
volatility clustering, with large returns expected to be followed by large returns 
and small returns by small returns (of either sign). This, together with the fact 
that autocorrelations of volatility measures decline very slowly [1] [2] [3], has the 
clear implication that long memory effects should somehow be represented in the 
process and this is not included in the geometric Brownian motion hypothesis. 
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One other hand, as pointed out by Engle[l], when the future is uncertain, 
investors are less hkely to invest. Therefore uncertainty (volatihty) would have 
to be changing over time. The conclusion is that a dynamical model for volatility 
is needed and a in Eq.(IT]), rather than being a constant, becomes a process 
by itself. This idea led to many deterministic and stochastic models for the 
volatility ([S] [5] and references therein). 

In a previous paper [7], using both a criteria of mathematical simplicity and 
consistency with market data, a stochastic volatility model was constructed, 
with volatility driven by fractional noise. It appears to be the minimal model 
consistent both with mathematical simplicity and the market data. This data- 
inspired model is different from the many stochastic volatility models that have 
been proposed in the literature. The model was used to compute the price 
return statistics and asymptotic behavior, which were compared with actual 
data. Deviations from the classical Black-Scholes result and a new option pricing 
formula were also obtained. The fractional volatility model, its predictions and 
comparison with data will be reviewed in Section 2. 

When this fractional volatility model was first presented, an interesting re- 
mark by an economist was " all right, the model seems to fit reasonably well the 
data, but where is the economics ? " . The same remark might be made about 
the simple geometric Brownian model, which does not even fit the data and is 
used by most of the of mathematical finance practitioners.. But, of course, our 
economist was right. The fractional volatility model seems to be a reasonable 
mathematical parametrization of the market behavior, but it is not sufficient 
to fit the data. One should also search for the mechanisms in the market that 
lead to the observed phenomena. No agent-based model can pretend to be the 
market itself, not even a realistic image of it. Nevertheless it may provide a 
surrogate model of the basic mechanics at work there. Therefore, the idea in 
this paper is to use stylized agent-based market models and find out which fea- 
tures of these models correspond to each one of the features of the mathematical 
parametrization of the data. 

2 The fractional volatility model. Induced volatil- 
ity, statistics of returns, option pricing and 
leverage 

The basic hypothesis for the model construction were: 

(HI) The log-price process log5t belongs to a probability product space 
17(8)57 of which the first one, fl, is the Wiener space and the second, O , 
is a probability space to be reconstructed from the data. Denote by a; G 
and uj G Q the elements (sample paths) in J7 and and by J^t and the 
(7— algebras in O and O generated by the processes up to t. Then, a particular 
realization of the log-price process is denoted 

log St (uj,uj ) 
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This first hypothesis is really not limitative. Even if none of the non-trivial 
stochastic features of the log-price were to be captured by Brownian motion, 
that would simply mean that St is a trivial function in Q. 

(H2) The second hypothesis is stronger, although natural. One assumes 

that, for each fixed lu , log St (^UjUJ^ is a square integrable random variable in 
ft. 

A mathematical consequence of hypothesis (H2) is that, for each fixed u , 

= fit[;^')dt + at(^»,Lo'^dB{t) (2) 

where /xt ^•,a; ^ and cTf ^•jW ^ are well-defined processes in fl. (Theorem 1.1.3 
in Ref.0) 

Recall that if {Xt,J-t} is a process such that 

dXt = fitdt + atdB (t) (3) 
with fit and at being J^f— adapted processes, then 



= \im^ {E{Xt+e-Xt)\Tt} 



— lim 



al[E{Xt+e-Xtf\Tt} 



(4) 



The process associated to the probability space could then be inferred 
from the data. According to Q, for each fixed cu realization in fl one has 

• ,Lo') = luni{£;(logS'i+e-log^t)'} (5) 

Each set of market data corresponds to a particular realization lu . Therefore, 
assuming the realization to be typical, the process may be reconstructed 
from the data by the use of (O. This data-reconstructed at process was called 
the induced volatility. 

For practical purposes we cannot strictly use Eq.(l5]) to reconstruct the in- 
duced volatility process, because when the time interval e is very small the 
empirical evaluation of the variance becomes unreliable. Instead, at was esti- 
mated from 

^.^= 1^ \ . var [log St) (6) 
\J-o - J-i\ 

with a time window jTo — Ti| sufficiently small to give a reasonably local char- 
acterization of the volatility, but also sufficiently large to allow for a reliable 
estimate of the local variance of log St- 

Once several data sets were analyzed[7], the next step towards obtaining a 
mathematical characterization of the induced volatility process was to look for 
scaling properties. It turned out that neither 

E\a{t + A)~a{t)\^A" (7) 
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nor 



E 



a{t + A)~a (t) 



.(0 - 

were good hypothesis for the induced volatihty process. It means that the 
induced volatihty process itself is not self-similar. 

If instead, one computes the empirical integrated log-volatility, one finds 
that it is well represented by a relation of the form 

t/s 

Y,^oga (nd) = (5t + R, it) (9) 

n=0 

the Re (t) process possessing very accurate self-similar properties. 

A nondegenerate process Xt, if it has finite variance, stationary increments 
and is self-similar 

Law (Xat) = Law {a" Xt) (10) 
must necessarily [9] have a covariance 

Cov {Xs,Xt) = \ + \t\'" -\s- E {Xl) (11) 

with < -ff < 1. The simplest process with these properties is a Gaussian 
process called fractional Brownian motion[TU]. with 

E [Bh m =0 it) Bh is)] = i + \s\'" - It - s\'"} (12) 

and, for i7 > ^ , a long range dependence 

oo 

J2CoviBHil),BHin + l)-B„in))=^ (13) 

n=l 

Therefore, mathematical simplicity suggested the identification of the R^r (t) 
process with fractional Brownian motion. 

Re it) =. kBn (t) (14) 

and, from the data, one obtains Hurst coefficients in the range 0.8 — 0.9. 
Finally one obtains the following fractional volatility model 

dSt = fiStdt + atStdB it) , 
log at = (3+l{BHit)-BHit-S)} 

fc is a volatility intensity parameter and d is the observation time scale. Notice 
that the volatility is not driven by fractional Brownian motion but by fractional 
noise, naturally introducing an observation scale dependence. 
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2.1 The statistics of price returns 

At each fixed time log at is a Gaussian random variable with mean (3 and vari- 
ance k^6^^~^. Then, 

/ ^ 1 /, X 1 f (logo--/9)M 

PS (a) = -PS (log a) = cxp -^p^^ (16) 



therefore 



with 



S' 



Ps (log ^) = dap, (a) (log ^) (17) 
' (log (t)-(^-^)(^ 



log —;- = — exp < 

^ StJ ^2na^ [T-t) ^ 



2a2 (T - t) 



(18) 

Thus, the effective probability distribution of the returns might depend both on 
the time lag A = T — f and on the observation time scale 5 used to construct 
the volatility process. That this latter dependence might actually be very weak, 
seems to be implied by comparison with the data from several markets. 

A closed-form expression for the returns distribution and its asymptotic 
behavior may be obtained, namely 

^'<^<^» = («-*<'"''-*'"r(^))|„, m 

with asymptotic behavior, for large returns 

Ps (r (A)) ~ -^=6-* ^ (20) 



with 



r(A) =log5T-log5t , = e^ A = r-i, A= ^''^^^ """^ 



2A6'2 



and 



Some illustrative comparisons with market data were performed. In Fig.l 
NYSE one-day data was used to fix the parameters of the volatility process. 
Then, using H = 0.83, k = 0.59, (3 = —5, (5 = 1, the one-day return distribution 
predicted by the model is compared with the data. The agreement is quite 
reasonable. For comparison a log-normal with the same mean and variance is 
also plotted in Fig.l. Then, in Fig. 2, using the same parameters, the same 
comparison is made for the A = 1 and A = 10 data. 
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Figure 1: One-day NYSE returns compared with the model predictions and the 
lognormal 

Fig. 3 shows a somewhat surprising result. Using the same parameters and 
just changing A from 1 (one day) to A = ^ (one minute), the prediction of the 
model is compared with one-minute data of USDollar-Euro market for a couple 
of months in 2001. The result is surprising, because one would not expect the 
volatility parametrization to carry over to such a different time scale and also 
because one is dealing with different markets. 

In Fig. 4 and Fig. 5 one sees the same one-day and one-minute return data 
discussed before, as well as the predictions of the model, both in semilogarithmic 
and loglog plots. 

As seen from Figs. 4 and 5, the exact result ^9]) or ((20|) resembles the 
double exponential distribution recognized by Silva, Prange and Yakovenko|llj 
as a new stylized fact in market data. The double exponential distribution has 
been shown, by Dragulescu and Yakovenko[T2, to follow from Heston's [13] 
stochastic volatility model. Notice however that our model is different from 
Heston's model in that volatility is driven by a process with memory (fractional 
noise). As a result, despite the qualitative similarity of behavior at intermediate 
return ranges, the analytic form of the distribution and the asymptotic behavior 
are different. 
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Figure 2: One and ten-days NYSE returns compared with the model predictions 
2.2 Option pricing 

New option pricing pricing formulas may be obtained from the model both in 
a simplified risk-neutral form or, more accurately, using fractional Malliavin 
calculus. Assuming risk neutrality [14j . the value V {St, at, t) of an option is the 
present value of the expected terminal value discounted at the risk-free rate 

V{St,at,t)^e-'^'^-''^ Jv{ST,aT,T)p{ST\St,at)dST (21) 

V {St, o'TjT) = max [0, S — K] and the conditional probability for the terminal 
price depends on St and at- K is the strike price, T the maturity time and St 
and at the price and volatility of the underlying security. 

In stochastic volatility models (with or without fractional noise) risk-neutrality 
is not an accurate assumption. Nevertheless it provides an approximate esti- 
mate of the deviations from Black-Scholes implied by the fractional volatility 
model. As in Hull and White [TS^, one uses the relation between conditional 
probabilities of related variables, namely 

p{ST\St,at) = p{ST\St,loga) p(loga\\ogat) d{loga) (22) 
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Figure 3: One-minute USD-Euro returns compared with the model predictions, 
with parameters obtained from one-day NYSE data 



log (T being the random variable 



T 



loga=^7— ^/ log a sds (23) 



that is, log cr is the mean volatility from time t to the maturity time T condi- 
tioned to an average value logct at time t. Finally the result for V {St,(7t,t) 

V{St,at,t) = St [aM {a,a,b) + bM {a,b,a)]-Ke~'''^^-^^ [aM {a, a, -b) - bM{a,-b, 

(24) 

-| poo poo 2 2 2 

M{a,a,b) = — / dy dxe-^e-'^(""+^) (25) 
27ra J_i Jq 

1 /2" . ^ «^ b 

= - — \ — / ax rsrfc — 

4a V ^ Jo aa; + I \ V2 V^x 

erfc is the complementary error function and a and b are 



(26) 
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Figure 4: Semilogarithmic and loglog plots of NYSE data 



In Fig. 6 one plots the option value surface for V {St, (Jt,t) in the range T—tG 
[5, 100] and S/K e [0.5, 1.5] as well as the difference {V {St,at,t) ~ C (5*, dt, i)) /K 
for A: = 1 and k = 2. The other parameters are fixed at ct = 0.01, r — 0.001, 6 = 
1, H = 0.8. To compare the predictions of the option pricing formula (j24p with 
the classical Black-Scholes (BS) result [1^ [17], the implied volatility required in 
BS to reproduce the same results was computed. This is plotted in the lower 
right panel of Fig. 6 which shows the implied volatility surface corresponding to 
V {St, (Jt, t) for k — 1. One sees that, when compared to BS, it predicts a smile 
effect with the smile increasing as maturity approaches. 



2.3 Leverage 

The following nonlinear correlation of the returns 



L{r) = {\r{t + T)\\{t))-{\r{t + T)\'){r{t)) 



(27) 



is called leverage and the leverage effect is the fact that, for t > Q, L (r) starts 
from a negative value and decays to zero whereas for r < it has almost 
negligible values. Fig. 7 shows L (r) computed for the NYSE index one-day data 
in the period 1966-2000. 
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Figure 5: Semilogarithmic and loglog plots of USD-Euro data 



The leverage behavior of the fractional volatility model will now be exam- 
ined. For this purpose it will be convenient to use the following integral repre- 
sentation of fractional Brownian motion [9] . 



Bh (t) = C 



H- 



dB{s)+ I {t- sf^^ dB{s) 





(28) 



Using this representation the fractional volatility model may be rewritten as 
dSt = fiStdt + atStdB(^Ht) 

log a* = f3 + k' f_^{t~sf-UB(^Hs) ^^^> 

where B^^-* (s) and B^^'' (s) different Brownian processes. In Fig. 8 one shows 
the leverage L (r) computed for the model ((29| with (3 and k chosen to match 
the statistical parameters of the NYSE index. Both the B'-^^ (s) ^ B^'^'> (s) 
and the _B'^' (s) = B^^'' (s) cases are considered. One sees that for S^^^ (s) ^ 
(s) there is no leverage effect, whereas for i?^^-' (s) = B^^'^ (s) an effect is 
found. Therefore one sees that, identifying the random generator of the log- 
price process with the stochastic integrator of the volatility, at least a part of 
the leverage effect is taken into account. 
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Figure 6: Option price and equivalent implied volatility in the "risk- neutral" 
approach to the stochastic volatility model 

3 Two agent-based studies 

Many factors play a role in a real market. To take into account all the factors in 
a model is neither possible nor, in many cases, illuminating. The objective is to 
isolate some of the more relevant mechanisms that presumably play a role in the 
market and, by stripping the model from other (inessential?) complications, to 
exhibit and understand the purified effect of these factors. As in other branches 
of science, the splitting apart of the dynamical components of a phenomena, 
may improve its understanding ^18j . 

Two stylized models will be considered. In the first the traders strategies 
play a determinant role. In the second the determinant effect is the limit-order 
book dynamics, the agents having a random nature. 

3.1 Agent strategies and market impact 

A market model with either random self-adapted strategies or fixed strategies 
was studied in detail in [19j . There it was found that the dominance of two 
types of strategies was to a large extent determined by the initial conditions. 
Different types of return statistics corresponded to the relative importance of 
either value investors or technical traders. The occurrence of market bubbles 
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Figure 7: Leverage for the NYSE one-day data in the period 1966-200 



also corresponded to situations where technical trader strategies were well rep- 
resented. 

Here, that model will be used for comparison purposes with the fractional 
volatility parametrization. The basic ingredients of the model are summarized 
below: 

One considers a set of investors playing against the market, that is, they 
have some effect on an existing market that is influenced by other factors (other 
investors and general economic effects). This assumption implies that in addi- 
tion to the impact of this group of investors on the market, the other factors 
are represented by a stochastic process. Therefore 

zt+i = f {zt,<^t) +Vt (30) 

represents the change in the log price (zf = log pt) with ujt being the total invest- 
ment made by the group of traders and r]t the stochastic process that represents 
all the other factors. No conservation law is assumed for the total amount of 
stock s and cash m detained by the group of traders. If pt is the price of the 
traded asset at time the purpose of the investors is to have an increase, as 
large as possible, of the total wealth mt+Pt x st at the expense of the rest of 
the market. 
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Figure 8: Leverage in the fractional volatility model 

The collective variable is z and each investor payoff at time t is 

A[') = (m« +p, X - (m« +P0 X (31) 
Market impact 

Let p be the price of a representative asset, z = log(p) and uit the total sum 
of the buying and selling orders (in money units) for the asset. Buying orders 
are positive and selling ones negative. The effect of these orders on the price 
change of the asset is called the market impact function. Let small orders have 
an impact according to the loglinear lawj^D] [H] 

zt+i ~ zt^'^+rjt (32) 

The constant A, is sometimes called the liquidity. Eg. ([5^ corresponds naturally 
to a first order expansion and satisfies the condition 

p(p(po, ^(1)), c.(2)) = ,,(1) + J^)) (33) 

which is expected to be valid for small orders. However, as pointed out by 
Zhang [22] there is experimental evidence that this is not an accurate represen- 
tation for large orders. Therefore a slightly different market impact function 
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will be used. The reasoning used to motivate it, has some relation to Zhang's 
although the result is somewhat different. 

When using Eq. (j32p in a discrete-time dynamical model we are somehow 
neglecting the fact that the market takes different times to fulfill (and to react 
to) small and large orders. Therefore this should be taken into account when 
reducing the dynamics to a sequence of equal time steps. In particular, the 
reaction of the market may be parametrized by a change in the A coefficient 
which, being also related to some stochastic process, may vary by a factor 
proportional to t". Taking the time t to fill an order to be proportional to its 
size, one obtains 

zt+i - zt = - — — 7a + r]t (34) 

This price impact function was first proposed in |19j with a = ^. For small 
orders it recovers the loglinear approximation and for very large orders (and 
a = i) Zhang's square root law. 

The agent strategies 

In first-order, two main types of informations are taken into account by the 
investors, namely the difference (misprice) between price and perceived value Vt 

^t- zt = \og{vt) ~ log(pt) (35) 

and the variation in time of the price (the price trend) 

zt - zt^i = log(pt) - log(pt_i) (36) 

Consider now a non-decreasing function f{x) such that /(— oo) = and /(oo) = 
1. Two useful examples are 

/i(^) = Oix) 

J2\^) — l+oxp(-/3£c) 

The information about misprice and price trend is coded on a four-component 
vector 7 

/ /(6 - zt)f{zt - Zt^l) \ 

fiCt ~ Zt) (1 - /(zt - Zt-l)) 

(i-/(6-^t))/(^t-^t-i) 
V (1 - /(6 - Zt)) (1 - !{zt - zt-i)) / 

The strategy of each investor is also a four-component vector a*^*) with entries 
— 1, 0, or 1. —1 means to sell, 1 means to buy and means to do nothing. 
Hence, at each time, the investment of agent i is a*-*-* • 7 . A fundamental (value- 
investing strategy) that buys when the price is smaller than the value and sells 
otherwise would be a^*^ = (1,1,— 1,-1) and a pure trend-following strategy 
would be a^*-* = (1,-1,1,-1) . In this setting the total number of possible 
strategies is B** = 81. The strategies will be labelled by numbers 
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(38) 



{^^ + 1) (39) 
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Therefore the fundamental strategy is strategy no. 72 and the pure trend- 
following one is no. 60. 

An evolution dynamics may be implemented in the model in the following 
way. After a number r of time steps, s agents copy the strategy of the s best 
performers and, at the same time, have some probability to mutate that strategy. 
This evolution aims at attaining the goal of improving gains, while at the same 
time allowing for some renewal of the strategies. The percentage of each strategy 
changes in time and one may find whether some of them become dominating or 
stable and when this may occur. 




Figure 9: Price, volatility and Ra (t) with equal amounts of fundamental and 
trend-following traders 

The model was run with different initial conditions and with or without 
evolution of the strategies. The following results were obtained in .19] : 

- When in the initial condition all traders have the fundamental strategy 
and evolution is activated, this strategy stays dominant, not being invaded by 
any other of the strategies that are created by the mutation process. There 
are however a few other strategies that, after being created, survive the selec- 
tion process. This is true for example for the strategies 45 = (0,1,-1,-1), 
18 = (-1,1,-1,-1), 73 = (1,1,-1,0) and 75 = (1,1,0,-1). These surviving 
strategies are however similar to the fundamental one. When there is domi- 
nance of the fundamental strategies, the price increments dp have a Gaussian 
distribution. 
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- The fundamental strategy ceases to be stable if it occurs in the initial 
condition in smaller amounts (< 40%). 

- For a completely random mixture of strategies in the initial condition, 
although the selection mechanism still favors at each evaluation cycle the best 
performers, the system never organizes itself to make the total payoff of this 
group of traders to grow, nor does a clear dominant strategy emerges. 

- When the simulation is run without evolution, with a fixed 50% of fun- 
damental strategies (no. 72) and 50% of trend-following ones (no. 60), one 
sees a large number of bubbles and crashes in the price evolution and the price 
increments distribution has fat tails. 

Because this last case is the one where the returns statistics is closer to the 
actual market data, it is here further analyzed to see whether it also displays 
the other features of the fractional volatility model. From typical simulation 
runs one computes 

t/s 

n=0 

and 

\R, {t + A)-R., {t)\ 

Fig. 9 shows a typical plot of the price process p{t), the volatility, {t) and 
¥^ {\Ra {t + A) — R„ {t)\} obtained from the model with equal amounts of funda- 
mental (1, 1, —1, —1) and trend-following (1,-1,1,-1) agents and no evolution. 
One notices the lack of scaling behavior of R^ {t) with an asymptotic exponent 
0.55, denoting the lack of memory of the volatility process. This might already 
be evident from the time behavior of R^ (t) in the lower left plot. Also, al- 
though the returns have fat tails in this case, they are of different shape from 
those observed in the market data. Similar conclusions are obtained with other 
combinations of agent strategies. In conclusion: It seems that the features of 
the fractional volatility model are not easily captured by a choice of strategies 
in an agent-based model. Notice however that what the fractional volatility 
model parametrizes is the bulk of the market data, that is, the behavior of the 
market in normal days. The agents reactions and strategies are very probably 
determinant during market crisis and market bubbles. 

3.2 A limit-order book dynamics model 

Here one considers a limit-order book where asks and bids arrive at random on 
a window [p — w,p + w] around the current price p. Every time a buy order 
arrives it is fulfilled by the closest non-empty ask slot, the new current price 
being determined by the value of the ask that fulfills it. If no ask exists when 
a buy order arrives it goes to a cumulative register to wait to be fulfilled. The 
symmetric process occurs when a sell order arrives, the new price being the bid 
that buys it. Because the window around the current price moves up and down, 
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asks and bids that are too far away from the current price are automatically 
eliminated. Sell and buy orders, asks and bids all arrive at random. The only 
parameters of the model are the width w of the limit-order book and the size n 
of the asks and bids, the sell and buy orders being normalized to one. 

The model was run for different widths w and liquidities n and, for compari- 
son with the fractional volatility model, one computes as before , X^nio ^ ("''^) ~ 
pt + R„ (t) and \R^ {t + A) - R„ {t)\. 



0.014 
0.012 





Figure 10: Price, volatility and R„ (t) in the limit order model 

Although the exact values of the statistical parameters depend on w and 
n, the statistical nature of the results seems to be essentially the same. Fig. 10 
shows typical plots of the price process p (t), the volatility, (t) and E {\Ra {t + A) — R„ {t)\} 
obtained for n — 2 and the limit-order book divided into 2w + 1 = 21 discrete 
price slots with Ap = 0.1. The scaling properties of R^ (t) are quite evident 
from the lower right plot in the figure, the Hurst coefficient being 0.96. Fig. 11 
shows the correlation and the pdf of the onc-timc returns. From these results 
one concludes that the main statistical properties of the market data (fast decay 
of the linear correlation of the returns, non-Gaussianity and volatility memory) 
are already generated by the dynamics of the limit-order book with random be- 
havior of the agents. This implies, as pointed out by some authors that in the 
past have considered limit order book models [23j [24j [25] [26] , that a large part 
of the market statistical properties (in normal business-as-usual days) depends 
more on the nature of the price fixing financial institutions than on particular 



17 




Figure 11: Linear correlation and pdf of returns in the limit order model 



investor strategies. 



4 Conclusions 

(a) The fractional volatility model provides a reasonable mathematical paramctriza- 
tion of the bulk market data, that is, it captures the behavior of the market in 
business- as-usual trading days. 

(b) A small modification of the original model, identifying the random gen- 
erator of the log-price process and the integrator of the volatility process, also 
describes, at least, a part of the leverage effect. 

(c) The market statistical behavior in normal days seems to be more influ- 
enced by the nature of the financial institutions (the double auction process) 
than by the traders strategies. Specific trader strategies and psychology should 
however play a role on market crisis and bubbles. 
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